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ABSTRACT 

We present a detailed analysis of the first very long baseline interferometry (VLBI) detection of 
the radio remnant of supernova 1987A. The VLBI data taken in 2007 and 2008 at 1.4 and 1.7 GHz, 
respectively, provide images sensitive to angular scales from O'.'l to 0'.'7, the highest resolution to date 
at radio frequencies. The results reveal two extended lobes with an overall morphology consistent with 
observations at lower resolutions. We find evidence of small-scale features in the radio shell, which 
possibly consist of compact clumps near the inner surface of the shell. These features have angular 
extent smaller than 0'.'2 and contribute less than 13% of the total remnant fiux density. No central 
source is detected in the VLBI images. We place a 3a flux density limit of 0.3 mJy on any pulsar or 
pulsar wind nebula at 1.7 GHz. 

Subject headings: Instrumentation: high angular resolution — supernovae: individual (SN 1987A) — 
ISM: supernova remnants 



1. INTRODUCTION 

Supernova (SN) 1987A in the Large Magellanic Cloud 
is the only nearby core-collapse supernova observed in 
the age of modern astrophysical instrumentation and 
has been intensively studied over the last two decades. 
Optical observations reveal a triple-ring nebula centered 
on the supernova, as part of the circumstellar material 
(GSM) expel led by the progenitor star in its late evo- 
lution stage (jCrotts fc HeathcQt3ll991f ). With a radius 
of 0'.'86, the nebula' s inner ring mark s the inner surface 
of the dense GSM (jPlait et al.lll995r ). Images from the 
Hubble Space Telescope (HST) s how hot spots em erging 
from the inner ring since 1995 (Law rence et al.l 12000). 
There are curre ntly about 30 optica l hot spots encircling 
the entire ring (jFrance et al.ll201Q[ ). which are believed 
to be clumps of inwar d-protruding GS M shocked by the 
supernova blast wave (jPun et al.ll2002f ). In recent years, 
the interaction between the blast wave and the inner ring 
has led to a rapid increase in the radio a nd soft X-ray 
emission (iZanardo et al.ll2010l : IPark et al.l[2007l) . provid- 
ing a good opportunity to study the shock interaction 
with the GSM. 

At radio frequencies, the initial outburst was de- 
tected with the M olonglo Observatory Synthesis Tele- 
scope (jTurtle et al.l 11987). the radio flux then decayed 
rapidly and re-emer ged in mid- 1990, markings; the birth 
of a radio remnant (jTurtle et al.lll99Qf ). Since then, the 
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source has been monitored regularly wit h the Australia 
Telescope Compact Array (ATCA) fsee iNg et all 120081 : 
IZanardo et al.l |2010L and references therein). Super- 
resolved images at 9 GHz reveal a two- lobe struc- 
ture of the radio shell (|Stavelev- Smith et al.l 119931 : 
iGaensler et al.l Il997f ) , which was later confirmed with 
diffract i on- limited imag:es at higher resolutions at 1 8 and 
36 GHz (Manchester et alj 120051 : iPotter et al.ll2009f ). 

High spatial resolution is the key to study the super- 
nova shock evolution and its interaction with the GSM. 
The first very long baseline interferometry (VLBI) ob- 
servations of SN 1987 A were made 5.2 days after the 
explosion. Although the radio emission was completely 
resolved and hence no detection was made, the obser- 
vations posted a lo wer limit of 1.9 x 10 ^kms~^ on the 
expansion velocity (i Jauncev et al.llT988f ). A subsequent 
attempt at VLBI in 2003 September also resulted in a 
non-detection (J. Lovell & R.W. Hunstead 2007, private 
communication). The first VLBI detection of the ra- 
dio remnant of SN 1987A was made with the Australian 
Long Baseline Array (L BAj^ in 2 007, and initial results 
have been reported by Ting av et al.l (120091) . In this Let- 
ter, we present a detailed analysis of the 2007 dataset, 
together with new observations in 2008, and present ev- 
idence of small-scale features in the radio shell. 

2. OBSERVATIONS AND DATA REDUCTION 

VLBI observations of SN 1987 A were carried out on 

2007 October 7 and 2008 November 26 at 1.4 and 
1.7 GHz, respectively, with about 10 hours integration 
time each, using part of the LBA including the ATCA 
(4 X 22 m antennas used as a tied array in 2007, and 5 
antennas in 2008), Parkes (64 m) and Mopra (22 m) tele- 
scopes of the Australia Telescope National Facility. The 

2008 epoch has additional data taken from the NASA 
DSS 43 antenna at Tidbinbilla (70 m) for three hours 
to boost the u-v coverage and sensitivity. The observa- 
tion parameters are listed in Table [TJ and the u-v cov- 

^ |http : //www . atnf . csiro . au/vlbi/overview/ 1 
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TABLE 1 

Observation and best-fit model parameters for 2007 and 
2008 VLBI OBSERVATIONS OF SN 1987A 



Parameter 



2007 



2008 



Observation Parameters 
2007 October 7 



Date 

Day since explosion 7531 
Center frequency (GHz) 1.382 

No. of stations 3 

Integration time (hr) 9 

Baseline (kA) 260-1250 

Baseline-hour 27 



2008 November 26 

7947 

1.666 

4 

10 

300-2510 
36 



Best-fit Truncated-shell Parameters^ 



Radius 

Half-opening Angle (°) 
Thickness (%) 



0.82 ±0.01 

23 ±4 
25 ±8 
0.81/25848 



0.84 ±0.01 
18.2 ± 1.2 
17.7 ± 1.4 
0.61/5122 



^ obtained with (i-^- weighted) fits; difference between the 
weighted and unweighted fits give the reported uncertainties. 

erage is plotted in Figure [H The 2007 and 2008 data 
are centered on 1.382 and 1.666 GHz, respectively, with 
four dual-polarization intermediate frequencies and a to- 
tal bandwidth of 64 MHz divided into four bands (except 
from the Tidbinbilla antenna, where only the left circu- 
lar polarization is available). A nearby calibrator, PKS 
0530—727, was observed every 10 minutes for phase ref- 
erencing purposes. 

The 2007 observations were made as part of an Express 
Production Real-time e-VLBI Service (EXPReS) demon- 
stration of a global real-time interferometer. The digital 
data for four dual-polarization 16 MHz bands were trans- 
ferred in real time from each antenna in Australia via a 
dedicated 1 Gbps lightpath to the correlation facility on 
the opposite side of the Earth (Joint Institute for VLBI 
in Europe [JIVE] in the Netherlands), using standard gi- 
gabit Ethernet network adapters and switches within the 
observatories. All four stokes parameters were correlated 
at the European VLBI Network (EVN) MarklV corre- 
lator at JIVE. Technical details on the li ghtpaths and 
correlation process have been reported bv iTingav et al.l 
(2009i ). The 2008 observations were performed in regular 
VLBI mode where the voltage data was recorded on disk, 
and subsequently correlated at the Curtin University fa- 
cility. 

All post-correlation calibration and processing were 
carried out with the AIPS package. Images were formed 
with natural weighting and then deconvolved using the 
CLEAN algorithm ( Clark 1 19801) . We performed self- 
calibration with a three-minute timescale. Since the cal- 
ibrator delivers high quality phase corrections, this of- 
fers only a minimal improvement to the images. Fig- 
ure[2K shows the final intensity maps of the 2007 and 2008 
observations, which have rms noise of 0.35 mJy beam~^ 
and 0.1 mJy beam" ^, and restoring beams of FWHM 
0'.'17 X 0'.'08 and 0'.'12 x 0'.'05, respectively. 

3. RESULTS AND DISCUSSION 

3.1. VLBI Results and Comparison with Observations 
at other Frequencies 

The VLBI images of SN 1987A in Figure [2^ clearly 
resolve the remnant structure and show a two-lobe mor- 
phology. The lobes are separated by 1'.'5 and are each 
extended with a size ~ 0'.'6 x 0'.'3, a few times larger 
the restoring beams. The eastern lobe has a higher flux 
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Fig. 1. u-v coverage of 2007 (black) and 2008 (red) VLBI ob- 
servations at 1.4 and 1.7 GHz, respectively, comparing to that of 
9 GHz ATCA observations (green). 

density than the western one by a factor of 2.5. There 
is some hint that the lobes consist of distinct peaks. A 
comparison between the 2007 and 2008 images show a 
similar overall morphology, but the lobes in the former 
appears to be more extended in R.A. and have slightly 
different orientations than those in the latter. There is 
also additional emission in 2007 to the south between the 
two lobes, which is not observed in 2008. We believe that 
these discrepancies are mostly due to artifacts in the im- 
age reconstru ction, as a result of s parse u-v sampling of 
the data (see lHevwood et "al]l2009f ). On the other hand, 
we argue that the two main lobes are physical, because 
they are present in both observations which have differ- 
ent hour angle coverage. The total cleaned flux density 
in the 2007 and 2008 images are 120 and 50mJy, re- 
spectively, much lower than the values 350 and 410 mJy 
obtained from ATCA monitori ng at a similar frequency 
in the same time period (Zanar do et al.|[2QTQl ). With a 
minimum baseline of ~300kA, the VLBI observations are 
insensitive to features larger than 0'.'7. Hence, our re- 
sult indicates that the majority of the remnant flux den- 
sity is contributed by smooth diffuse emission of at least 
a comparable scale, implying a substantial thickness of 
the radio shell. 

In Figure [213 the 2008 VLBI image is compared to lower 
resolution ATCA images at 9 and 36 GHz t aken at sim- 
ilar e pochs (2008 October) (jNg et al.ll207)8l : iPotter et all 
l2009l ). Since all these observations are phase-referenced 
to the same nearby calibrator, we estimate that the as- 
trometric accuracy should be ~ 10 mas and therefore did 
not perform any manual alignment. The lobes found in 
the VLBI images closely match the brightest part of the 
shell at 9 and 36 GHz. However, regions of lower bright- 
ness temperature tracing the remainder of the remnant, 
seen in the latter two images, are not detected in the 
VLBI observations. This is likely because the large-scale 
diffuse emission is resolved out in the VLBI images, and 
any compact features (see §3.31 below) fall below the de- 
tection sensitivity. Figure [JJd also presents a comparison 
of the VLBI result to Ha and X-ray images obtained 
from the B.ST and the Chandra X-ray Observatory on 
2007 May and 2008 April, respectively (J.C.S. Pun, pri- 



VLBI Detection of SN 1987A 



3 




5^35™28'2 28"! 28'0 27^9 27^8 5'^35"'28'2 28'l 28'0 27^9 27^8 

RA (J2000) RA (J2000) 



(b) 



-69 16 12.5 - 




28'1 28'0 27'9 
RA (J2000) 



28'1 28"0 27'9 
RA (J2000) 



Fig. 2. (a) VLBI images of SN 1987A at 1.4 and 1.7 GHz, taken on 2007 October 7 and 2008 November 26, respectively. The contours 
correspond to 1, 4, SmJybeam"-*^ in 2007, and 0.5, 1.5, 3 and 5mJybeam~-'^ in 2008. The restoring beams are shown in lower left, and 
the scale bars are in units of Jybeam"-*^. (b) Comparison of the 2008 VLBI results to: super-resolved 9 GHz ATCA image taken on 2008 
October 11 (Ng et al. ■200& V diffraction-hmited 36 GHz ATCA image taken on 2008 October 7-12 (Potter et al. 2009), HST Ho; images 
taken on 2007 May 11 and 2000 November 14 (inset), and deconvolved Chandra X-ray image in 0.3 — lOkeV band, taken on 2008 April 
28-29 (Ng et al. 20091). 

The model parameters include the center position, flux, 
mean shell radius (averaged between the inner and outer 
shells), half-opening angle, shell thickness, and a linear 
gradient in the surface emissivity th at accounts for the 
observed asymmetry (see Figure 3 in iNg et alJ l2QQ8l for 
a detailed definition of all these parameters,). Given the 
limited u-v coverage of the da ta, we fixed the gradient 
using the 9 GHz results (Potte r et al.l [20 09) in order to 
improve the fit stability. The level of gradient is fixed 
at 36% (ratio between the surface brightness at the edge 
and at the center) with a position angle at 109° (mea- 
sured from the observer's line-of-sight projected on the 
shell's equatorial plane). All other parameters, including 
position, fiux, radius, half-opening angle and thickness, 
are free to vary in the fit. The values are calculated 
by weighting the fit residuals with the measurement un- 



vate communication: INg et al.l[2009h . These images were 
fitted to a simple torus to determine the remnant center, 
then aligned with t he VLBI image usin g the supernova 
position reported by iPotter et al.r(|2009l ) . There appears 
to be no detailed correlation of the brightness distribu- 
tion at small scales in the radio, optical and X-ray bands. 
Moreover, the images indicate that the radio shell has a 
comparable size to the optical inner ring and the X-ray 
shell. Th e latte r result is consistent with the findings by 
INg et al.l(|20Q9h . 

3.2. Fourier Modeling 

To quantify the remnant geometry, we employ ed the 
Fourier modeling technique described bv iNg et"al] ([2008) 
to fit the VLBI visibility data with a truncated-shell 
model tilted at 43?4 to the line-of-sight (^Pun .20071 
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Fig. 3. — Images made from simulations of the 2008 VLBI visibility data, based on the best-fit truncated-shell models on the 2008 VLBI 
data (left) and on the 9 GHz ATCA data (right). The u-v sampling of the simulations are identical to the 2008 VLBI observation. The 
images are noise- free, and the contours are from the 2008 VLBI image in Figure [2] 



certainties of each data point, which are estimated from 
the scatter of the visibihty amphtudes among the four 
channels, because information on individual measure- 
ment errors were not available. Table [1] summarizes the 
main results: the best-fit models have radii 0^82 and 
0'/84 in 2007 and 2008, respectively, half-opening angle 
~ 20° and thickness ~ 20% of the radius. We note 
that the measured thickness should be considered only 
as a lower limit, since the data are insensitive to large- 
scale structure. The small reduced values suggest that 
the measurement uncertainties are likely overestimated. 
To determine the confidence intervals of the best-fit pa- 
rameters, we scaled the uncertainties such that the re- 
duced = 1, and found very small formal statistical 
errors. This implies that the systematic errors domi- 
nate, which could be introduced by calibration errors or 
by the weighting scheme of the fits. While the former is 
difficult to estimate, we attempted to quantify the lat- 
ter with simple least-square fits (i.e. without weighting). 
The difference in the model parameters give the uncer- 
tainties reported in Table [H 

Based on the best-fit model parameters for 2008, we 
simulated a visibility data set using the u-v sampling 
of the 2008 VLBI data. The simulated data were then 
imaged and deconvolved identically as the real observa- 
tions, and the resulting image is shown in Figure [3^. The 
truncated-shell model successfully captures the overall 
shell morphology, but fails to reproduce the double-peak 
feature of the eastern lobe. A direct comparison between 
the best-fit parameters in 2007 and 2008 indicates no sig- 
nificant change in the shell's geometry. The difference in 
radius gives an expansion veloci ty of 4000 ±3000 km s~^, 
consistent with the result from iNg et al.l (^2008), within 
the sizeable error due to the relatively short time pe- 
riod between the observations. The VLBI images show 
no obvious point source at the remnant center, and the 
2008 VLBI image gives a 3cr upper limit of 0.3 mJy on 
the flux density of any possible central object, such as a 
radio pulsar or pulsar wind nebula at 1.7 GHz. 

As listed in Table [TJ the best-fit radii are significantly 
smaller than the value 0'.'89 obt ained fro m the ATCA 
9 GHz observations in 2008 (Potte r et all [2009). To de- 
termine if the discrepancy could be caused by the lim- 
ited u-v coverage of the VLBI data, we simulated an im- 
age with the best-fit truncated-shell model of the 9 GHz 



data, and using the 2008 VLBI observations' u-v cover- 
age. The result is shown in Figure [Sb, clearly indicating 
a larger shell size than observed. This suggests that if 
the 9 GHz model were the correct description of the rem- 
nant geometry, then the VLBI observations should be 
able to recover the geometric parameters. We compare 
in Figure the mean shell radius obtained from Fo urier 
modeling of different radio observations (Ng et al. 2008|; 
iPotter et al.l l2009l ). and found that the best-fit radius 
appears to be largest at 9 GHz and smallest at 1.4 GHz, 
while the sizes at 18 and 36 GHz are in between. 

Discrepancy in shell radius at different frequencies has 
been claimed for the radio remnant of SN 1993 J, with 
a progressively sm aller radius a t 5 GHz than at 1.6 GHz 
since day 1500 (M arcaide et al.li2009 )Fl A few possible 
explanations have been proposed, including difference in 
electron synchrotron lifetimes, a radially increasing mag- 
netic field, or freq uency- dependen t changes in the opac- 
ity of the medium ("Marcaid e et al.l l2009) . For our case of 
SN 1987A, the non-monotonic change of the mean shell 
radius with frequency presents a more complex picture, 
suggesting that the above physical explanations may not 
apply. In Figure |4)3 we plot the deviation of radius from 
a linear fit to the 9 GHz radius after day 4000 against 
the theoretical (i.e. diffraction-limited) angular resolu- 
tion. We found a trend of decreasing radius as the beam 
size gets smaller, independent of the observing frequency. 
We have examined whether biases in our fitting technique 
may give rise to such a result, but as the fits are carried 
out in the Fourier domain rather than the image domain, 
our procedure is relatively immune to this. Therefore, 
we suggest that measurements at higher spatial resolu- 
tions are biased by unmodeled components in the radio 
emission, likely to be bright knots near the shell's inner 
surface. 

To test this idea, we tried applying the Fourier model- 
ing to the 18 and 36 GHz ATCA data, but restricted on 
the u-v range < 180kA, corresponding to the same cov- 
erage as the 9 GHz ATCA observations. In all cases, this 
results in a ~ 5% increase in the best-fit radii when com- 
pared to fits on the full dataset, giving values closer to the 
9 GHz results, thus, providing further support to the pic- 



However, we note that IBartel et al.l (|2002l ) did not find such 
a correlation using the same dataset. 
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Fig. 4. — (a) Radii of supernova remnants 1987A and 1993 J obtained from different radio observations (Ng et all 120081 : IPotter et al.l 
[20091 : [Marcalde et al. 2009), assuming distances of 51.4kpc and 3.96 Mpc, respectively (Bartel et al. 2007). (b) Deviation of SN 1987A's 
shell radius versus the theoretical resolution of different observations. The deviations are estimated from a linear fit to the 9 GHz radius 
beyond day 4000, which is shown by the solid line in (a). 

ture above. As a final check, we simulated visibility data 
based on a toy model of a truncated shell plus few point 
sources. The shell has a radius of 0'.'9 and 10% thickness, 
and we added four point sources, two at each side, at a 
distance O'.'S from the center and with a total flux den- 
sity of 15% of the remnant. Based on this toy model, we 
simulated two visibility data sets using the u-v sampling 
of the 2008 VLBI and ATCA 9 GHz observations, then 
applied the Fourier modeling as in the real observations. 
This gives radii of 0'.'85 and O'.'SS for the two simulated 
data sets, respectively, confirming that our radius mea- 
surements made with a truncated shell could be biased 
by compact features in the remnant, especially at high 
spatial resolution . This result is consistent with previ- 
ous findings (Man chester et al.l [2002: Ng et al. 2008), in 
which the radius is always smaller when fitted with a sin- 
gle shell than with a shell plus point sources. We note 
that although previous results at 9 GHz based on Fourier 
modeling (e.g. 'Ng et al.|[2QQ8h are likely biased too, our 
simulations suggest that this is at a much lower level than 
the VLBI case. Also, the bias is systematic, therefore, 
it should have a minimal effect on the expansion rate 
measurements. 



3.3. Small-scale Features 

Results from the Fourier modeling above present ev- 
idence of small-scale features in the radio remnant of 
SN 1987A. As hinted by the double-peak structure of the 
radio lobes seen in Figure [2^, these features likely con- 
sist of discrete clumps near the inner surface of the shell. 
To constrain their flux density and spatial extent, we as- 
sume an extreme case that all diffuse emission is resolved 
in the VLBI observations, such that the structure seen 
in the images represents only the compact emission. We 
fitted the 2008 VLBI visibility data with a simple model 
of four Gaussian blobs, two at each lobe, and found that 
the clumps have size smaller than 0'.'2, with a distance 
0'.'78-0'.'83 from the remnant center, and a total flux den- 
sity < 50 mJy. The size limit is only marginally larger 
than the restoring beam, and converts to a physical ex- 
tent < 1.5 X 10^'^ cm 0.05 pc) at the distance 51.4 kpc. 
This is comparable to the scale < 10^^ cm of the optical 



hot spots, for which hydrodynamic simulations suggest 
a physical siz e of 2 x 10^^ cni , too small to be resolved 
by the HST (iPun et al.l 120021 ). While it is tempting to 
associate the radio clumps with the optical hot spots, 
which are dense gas protruding inward from the inner 
ring. Figure [SJd shows no correlation between the radio 
and optical brightness distribution. Also, the first optical 
hot spot emerged in the nort heast with a position angle 
29° (see inset of Figure [SJd; Lawre nce et al.|[2QQQf ). sig- 
nificantly offset from the radio peak. Alternatively, the 
compact radio emission could indicate sites of enhanced 
cosmic ray acceleration, either with magnetic field am- 
plification or with a local field orientation parallel to the 
shock normal. Hence, tracing the evolution of the clumps 
in future observations may give direct constraints on the 
acce leration time scale an d the magnetic field strength 
(see lUchivama et al.l I2007D . 

4. CONCLUSION 

In this paper, we present a detailed study of VLBI ob- 
servations of SN 1987A taken in 2007 and 2008 at 1.4 
and 1.7 GHz, respectively. This gives the first VLBI de- 
tection of the radio remnant, revealing extended lobes 
that are consistent with the overall shell morphology in 
previous radio studies. We found that the Fourier mod- 
eling made with a simple truncated shell gives a smaller 
mean shell radius for the VLBI data than for other radio 
observations at lower spatial resolutions. This suggest 
small-scale features in the radio emission that are not 
captured by the smooth model. These features are likely 
to be discrete clumps near the shell's inner surface. We 
obtain a limit of 0'.'2 on their spatial scales, with a flux 
density limit 50mJy, corresponding to ~ 13% of that of 
the remnant. 

Although only part of the shell is detected in our VLBI 
observations, this study demonstrated VLBI imaging as 
a powerful tool to trace the fine structure of the radio 
emission. Further observations utilizing the full capa- 
bility of the LBA, including more antennas and wider 
bandwidth, will offer images with better sensitivity and 
resolution. In addition, the on-going interaction between 
the supernova blast wave and the dense GSM has induced 
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a rapid brightening of the remnant. Altogether, we ex- 
pect future VLBI images to reveal more clumps around 
the shell, and potentially their substructure. The new 
observations will also provide a direct measurement of 
the remnant's expansion rate, tracing the shock evolu- 
tion in the CSM. Finally, we note that if in future the 
clumps start to dominate the remnant flux, this will have 
a significant impact on the Fourier modeling of the radio 
data. In particular, the east-west asymmetry of the shell 
could depend sensitively on the brightness of individual 
clumps, as a result, the simple linear gradient approxi- 
mation used by iNg et al.l ()2QQ8f ) may no longer hold. 
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